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Abstract 



We report a measurement of the W boson mass based on an integrated lumi- 
nosity of 82 pb _1 from pp collisions at yfs = 1.8 TeV recorded in 1994-1995 
by the D0 detector at the Fermilab Tevatron. We identify W bosons by their 
decays to ev and extract the mass by fitting the transverse mass spectrum 
from 28,323 W boson candidates. A sample of 3,563 dielectron events, mostly 
due to Z — > ee decays, constrains models of W boson production and the de- 
tector. We measure M w = 80.44±0.10(stat)±0.07(syst) GeV. By combining 
this measurement with our result from the 1992-1993 data set, we obtain 
M w = 80.43 ±0.11 GeV. 
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In the standard model of the electroweak interactions (SM) [T|], the mass of the W boson 
is predicted to be 



In the "on-shell" scheme @ cos 6 W = Mw/Mz, where Mz is the Z boson mass. A measure- 
ment of M w , together with M z , the Fermi constant (Gp), and the electromagnetic coupling 
constant (a), determines the electroweak radiative corrections Ar experimentally. Purely 
electromagnetic corrections are absorbed into the value of a by evaluating it at Q 2 = M§. 
The dominant contributions to Ar arise from loop diagrams that involve the top quark and 
the Higgs boson. If additional particles which couple to the W boson exist, they will give 
rise to additional contributions to Ar. Therefore, a measurement of My/ is one of the most 
stringent experimental tests of SM predictions. Deviations from the predictions may indi- 
cate the existence of new physics. Within the SM, measurements of Myy and the mass of 
the top quark constrain the mass of the Higgs boson. 

This Letter reports a precise new measurement of the W boson mass based on an in- 
tegrated luminosity of 82 pb -1 from pp collisions at y/s = 1.8 TeV, recorded by the D0 
detector || during the 1994-1995 run of the Fermilab Tevatron. A more complete account 
of this analysis can be found in Refs. Previously published measurements |7|-[T3|, when 

combined, determine the W boson mass to a precision of 125 MeV. 

At the Tevatron, W bosons are produced mainly through qq annihilation. We detect 
them by their decays into electron-neutrino pairs, characterized by an isolated electron |TJ| 
with large transverse momentum (px) and significant transverse momentum imbalance ($t)- 
The $t is due to the neutrino which escapes detection. Many other particles of lower 
momenta, which recoil against the W boson, are produced in the breakup of the proton and 
antiproton. We refer to them collectively as the underlying event. 

At the trigger level we require ^ > 15 GeV and an energy cluster in the electromag- 
netic (EM) calorimeter with pt > 20 GeV. The cluster must be isolated and have a shape 
consistent with that of an electron shower. 

During event reconstruction, electrons are identified as energy clusters in the EM 
calorimeter which satisfy isolation and shower shape cuts and have a drift chamber track 
pointing towards the cluster centroid. We determine their energies by adding the energy 
depositions in the first ~ 40 radiation lengths of the calorimeter in a window, spanning 0.5 
in azimuth ((f)) by 0.5 in pseudorapidity (rj) UTS] , centered on the highest energy deposit 
in the cluster. Fiducial cuts reject electron candidates near calorimeter module edges and 
ensure a uniform calorimeter response for the selected electrons. The electron momentum 
(p(e)) is determined by combining its energy with its direction which is obtained from the 
shower centroid position and the drift chamber track. The trajectories of the electron and 
the proton beam define the position of the event vertex. 

We measure the sum of the transverse momenta of all the particles recoiling against the 
W boson, ut = J2i Ei sin 8iUi, where E± is the energy deposition in the i th calorimeter cell 
and 9i is the angle defined by the cell center, the event vertex, and the proton beam. The 
unit vector Wj points perpendicularly from the beam to the cell center. The calculation of «y 
excludes the cells occupied by the electron. The sum of the momentum components along 




(1) 
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the beam is not well measured because of particles escaping through the beam pipe. From 
momentum conservation we infer the transverse neutrino momentum, Pt{v) — — Pr( e ) -«Tj 
and the transverse momentum of the W boson, pr(W) = —ut- 

We select a W boson sample of 28,323 events by requiring priy) > 25 GeV, ut < 15 GeV, 
and an electron candidate with \rj\ < 1.0 and pr(e) > 25 GeV. 

Since we do not measure the longitudinal momentum components of the neutrinos 
from W boson decays, we cannot reconstruct the ev invariant mass. Instead, we ex- 
tract the W boson mass from the spectra of the electron p? and the transverse mass, 



rriT = y 2pr{e)pT{y){\ — cos A0), where A0 is the azimuthal separation between the two 
leptons. We perform a maximum likelihood fit to the data using probability density func- 
tions from a Monte Carlo program. Since neither ttit nor pr(e) are Lorentz invariants, we 
have to model the production dynamics of W bosons to correctly predict the spectra. The 
nix spectrum is insensitive to transverse boosts at leading order in Pt{W) / 'Myy and is there- 
fore less sensitive to the W boson production model than the pr(e) spectrum. On the other 
hand, the spectrum depends strongly on the detector response to the underlying event 
and is therefore more sensitive to detector effects than the py(e) spectrum. 

Z bosons decaying to electrons provide an important control sample. We use them to 
calibrate the detector response to the underlying event and to the electrons, and to constrain 
the model for intermediate vector boson production used in the Monte Carlo simulations. 

AZ->ee event is characterized by two isolated high-pr electrons. We trigger on events 
with at least two EM clusters with pt > 20 GeV. We define two samples of Z — ► ee decays in 
this analysis. For both Z samples, we require two electron candidates with px > 25 GeV. For 
sample I, we loosen the pseudorapidity cut for one of the electrons to \rj\ < 2.5. This selection 
accepts 2,341 events. For sample II, we require both electrons with \rj\ < 1.0 but allow 
one electron without a matching drift chamber track. Relaxing the track requirement for 
electrons with \rj\ < 1.0 increases the efficiency without a significant increase in background. 
Sample II contains 2,179 events of which 1,225 are in common with sample I. 

For this measurement we developed a fast Monte Carlo program that generates W and 
Z bosons with the rapidity and pt spectra given by a calculation using soft gluon resumma- 
tion TjJ and the MRSA' ]T7|] parton distribution functions. The line shape is a relativistic 
Breit-Wigner, skewed by the mass dependence of the parton luminosity. The measured 
intrinsic widths [|H|,|n| are used. The angular distribution of the decay electrons includes 



a pr(W)- dependent 0(a s ) correction |p0| . The program also generates W — > ewy [|21 
Z — > ee7 f21| , and W — > tu — > evvv decays. 



The program smears the generated p(e) and ut vectors using a parameterized detec- 
tor response model and applies inefficiencies introduced by the trigger and event selection 
requirements. The model parameters are adjusted to match the data and are discussed 
below. 

The energy resolution for electrons with \rj\ < 1.0 is described by sampling, noise, and 



constant terms. In the Monte Carlo simulation we use a sampling term of 13%/y7>r/GeV, 
derived from beam tests. The noise term is determined by pedestal distributions derived 
from the W data sample. We constrain the constant term to cem = l-15^Qjg% by requiring 
that the width of the dielectron invariant mass spectrum predicted by the Monte Carlo 
simulation is consistent with the Z data (Fig. [1]). 



6 



Beam tests show that the electron energy response of the calorimeter can be parameter- 
ized by a scale factor «em and an offset Sem- We determine these in situ using 7r° — > 77, 
J/ip -> ee, and Z -> ee decays. We obtain <5 E m = — 0.16±g;^ GeV and a EM = 0.9533±0.0008 
by fitting the observed mass spectra while constraining the resonance masses to their mea- 



sured values fO|pS| . The uncertainty on c^em is dominated by the finite size of the Z sample. 
Figure [l| shows the observed dielectron mass spectrum from sample II, and the line shape 
predicted by the Monte Carlo simulation for the fitted values of cem, oem, and <5em- 

We calibrate the response of the detector to the underlying event, relative to the EM 
response, using sample I. The looser rapidity cut on one electron brings the rapidity dis- 
tribution of the Z bosons closer to that of the W bosons, since there is no rapidity cut on 
the unobserved neutrino in W events. In Z — > ee decays, momentum conservation requires 
pr{ee) = —ut, where priee) is the sum of the two electron px vectors. To minimize sensi- 
tivity to the electron energy resolution, we project ut and pr(ee) on the inner bisector of 
the two electron directions, called the 77-axis (Fig. Q). We call the projections u v and p n {ee). 

Detector simulations based on the geant program p3j predict a detector response to the 
recoil particle momentum of the form R rec = a rec + /5 rec ln(p T /GeV). We constrain a rec and 
Prec by comparing the mean value of u v + p v (ee) with Monte Carlo predictions for different 
values of the parameters. We measure a rec = 0.693 ± 0.060 and j3 rec = 0.040 ± 0.021 with a 
correlation coefficient of —0.98. 

The recoil momentum resolution has two components. We smear the magnitude of 
the recoil momentum with a resolution of s rcc / y^r/GeV. We describe the detector noise 
and pile-up, which are independent of the boson p? and azimuthally symmetric, by adding 
the p't from a random pp interaction, scaled by a factor a m b, to the smeared boson p?. 
To model the luminosity dependence of this resolution component correctly, the sample of 
pp interactions was chosen to have the same luminosity spectrum as the W sample. We 
constrain the parameters by comparing the observed rms of w^/Rrec + Pn( ee ) with Monte 
Carlo predictions and measure s rec = 0.49 ±0.14 and a m b = 1.032 ±0.028 with a correlation 
coefficient of —0.60. Figure |2| shows a plot of u v /R rec + p v (ee). 

Excluding the cells occupied by the electrons, the average transverse energy flow, St = 
J2i Ei sin 9i, is 7.7 GeV higher for the W sample than for the Z sample. This bias is caused 
by requiring the identification of two electrons in the Z sample versus one in the W sample. 
The larger energy flow translates into a slightly broader recoil momentum resolution in the 
W sample. We correct a m b by a factor 1.03 ± 0.01 to account for this effect in the W boson 
model. 

Backgrounds in the W sample are W — ^ tu — >• evvv decays (1.6%), hadrons misidentified 
as electrons (1.3%±0.2%), Z -> ee decays (0.42%±0.08%), and W -> rv -> hadrons ± X 
decays (0.24%). Their shapes are included in the probability density functions used in the 
fits. 

The fit to the m T distribution (Fig. |3|(a)) yields M w = 80.44 GeV with a statistical 
uncertainty of 70 MeV. A Kolmogorov-Smirnov (KS) test gives a confidence level of 28% 
that the parent distribution of the data is the probability density function given by the Monte 
Carlo program. A x 2 test gives \ 2 = 79.5 for 60 bins which corresponds to a confidence 
level of 3%. The fit to the p T (e) distribution (Fig. |(b)) yields M w = 80.48 GeV with a 
statistical uncertainty of 87 MeV. The confidence level of the KS test is 83% and that of the 
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X 2 test is 35%. 

We estimate systematic uncertainties on Mw from the Monte Carlo parameters by vary- 
ing the parameters within their uncertainties. Table I summarizes the uncertainties in the 
W boson mass. In addition to the parameters described above, the calibration of the elec- 
tron polar angle measurement contributes a significant uncertainty. We use muons from pp 
collisions and cosmic rays to calibrate the drift chamber measurements and Z —>■ ee decays 
to align the calorimeter with the drift chambers. Smaller uncertainties are due to the re- 
moval of the cells occupied by the electron from the computation of ut, the uniformity of 
the calorimeter response, and the modeling of trigger and selection biases ||. 

The uncertainty due to the model for W boson production and decay consists of sev- 
eral components (Table I). We assign an uncertainty that characterizes the range of varia- 
tions in Mw obtained when employing several recent parton distribution functions: MRSA', 
MRSD— ' |1, CTEQ2M H, and CTEQ3M @. We allow the pr(W) spectrum to vary 



within constraints derived from the />r(ee) spectrum of the Z data |J and from Aqcd |22 



The uncertainty due to radiative decays contains an estimate of the effect of neglecting 
double photon emission in the Monte Carlo simulation . 

The fit to the tjit spectrum results in a W boson mass of 80.44 ± O.lO(stat) ± 
0.07(syst) GeV and the fit tothep T (e) spectrum results in 80.48±0.11(stat)±0.09(syst) GeV. 
The good agreement of the two fits shows that our simulation models the W boson produc- 
tion dynamics and the detector response well. We have performed additional consistency 
checks. A fit to the p T (v) distribution yields M w = 80.37 ± 0.12(stat) ± 0.13(syst) GeV, 
consistent with the tut and pr(e) fits. Fits to the data in bins of luminosity, 0(e), 77(e), and 
ut do not show evidence for any systematic biases. 

We combine the results from the tut fit and the data collected by D0 in 1992-1993 |J to 
obtain M\y = 80.43 ±0.11 GeV. This is the most precise measurement of the W boson mass 
to date. This result is in agreement with the prediction of 80.278 ±0.049 GeV from a global 
fit to electroweak data \T§\. Using Eq. [I] we find Ar = —0.0288 ± 0.0070, which establishes 



the existence of electroweak corrections to Mw at the level of four standard deviations. 
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TABLE I. Uncertainties in the W boson mass measurement in MeV, rounded to the nearest 5 



MeV. 



Source of uncertainty 


rriT fit 


pr(e) fit 


W sample size 


70 


85 


Z sample size (ckem) 


65 


65 


Total statistical uncertainty 


95 


105 


calorimeter linearity (<5em) 


20 


20 


calorimeter uniformity 


10 


10 


electron resolution (cem) 


25 


15 


electron angle calibration 


30 


30 


electron removal 


15 


15 


selection bias 


5 


10 


recoil resolution (a m b,Srec) 


25 


10 


recoil response (a rcc ,(3 vcc ) 


20 


15 


Total detector systematics 


60 


50 


Backgrounds 


10 


20 


Pt(W) spectrum 


10 


50 


parton distribution functions 


20 


50 


parton luminosity 


10 


10 


radiative decays 


15 


15 


W boson width 


10 


10 


Total W boson production and decay model 


30 


75 


Total 


115 


140 
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FIG. 1. The dielectron invariant mass distribution of the Z data for sample II (•). The solid 
line shows the fitted signal plus background shape and the small shaded area the background. The 
arrows indicate the fit window. 




FIG. 2. The definition of the 77-axis (left). 
(•) and simulation ( — ). 



The plot of u^/Rrec + Pr/(ee) (right) for the data 
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FIG. 3. Spectra of (a) niT and (b) pr(e) f 
(shaded). The arrows indicate the fit windows 




p T (e) (GeV) 



the data (•), the fit ( — ), and the backgrounds 



